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Neutron stars are discussed as laboratories of physics of strong gravi- 
tational fields. The mass of a neutron star is split into matter energy and 
gravitational field energy contributions. The energy of the gravitational 
field of neutron stars is calculated with three different approaches which 
give the same result. It is found that up to one half of the gravitational 
mass of maximum mass neutron stars is comprised by the gravitational 
field energy. Results are shown for a number of realistic equations of state 
of neutron star matter. 

PACS numbers: 04.40.Dg, 97.60.Jd 

1. Introduction 

Gravity makes all selfgravitating fluid bodies spherical. Thus the spher- 
ical symmetry of non-rotating neutron stars, which are real objects, is their 
intrinsic property, not an idealization. Neutron stars are relativistic objects, 
with mass of order IMq enclosed within 10-20 km. Gravity of such a com- 
pact mass distribution is strong enough to require relativistic description. 
In the following we focus on gravitational field contribution to the neutron 
star masses in General Relativity. 

The problem of the energy distribution of gravitational field remains 
an unsettled issue in general relativity mainly because of the lack of 
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proper energy-momentum tensor for gravitational field. However, for static, 
spherically symmetric space-time no controversy exists and relevant matter 
and gravity energies can be unambiguousely denned [HE]- Our main concern 
here is to study if neutron stars could shed some new light on the problem 
of gravitational field energy and its localization. 

To calculate both matter and gravitational field contributions to the 
neutron star mass we solve the Einstein's field equations for the static, 
spherically symmetric metric and a given equation of state of dense matter. 
The metric is 

ds 2 = e u c 2 dt 2 - e x dr 2 - r 2 (d6 2 + sin 9 2 d(f> 2 ) , (1) 

where v, A are functions of the radial coordinate r. 

Dense matter in neutron stars is thought to be in a liquid state, hence 
it is described by the perfect fluid energy-momentum tensor, 

T"" = (pc 2 + P)u»u v - Pg^, 

where pc 2 is the energy denisty, P is the pressure and is the four velocity. 
For the neutron star matter both energy density and pressure are uniquely 
determined by the baryon number density, p = p(hb),P = P{ub), hence 
the equation of state is of barotropic form, P = P{p). In this paper we 
employ a number of realistic equations of state, described in some details in 
Sect. E| 

Using the Einstein's equations one finds that functions A and v satisfy 
ordinary differential equations 



x// x 1 / a 8 vrGe A r 2 p(r) \ . . 

A'(r) = - \l - e x + ^^ij (2a) 

„ , 1 / \ 8TiGe x r 2 P(r)\ 
u'(r) = - f -1 + e A + - 4 (2b) 

p '(r) = ~(P(r) + pc 2 y . (2c) 

The metric (1) has to be continuous in the whole space. For this to 
happen, at the stellar boundary, r = R, the internal metric has to go over 
into the external Schwarzschild metric. There is a vacuum outside the star 
and the matter energy-momentum tensor vanishes. The metric functions at 
the boundary should satisfy the condition e u ^ = e~ x ^ = 1 — 2GM/Rc 2 , 
where M is the mass of the star. 
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From Eqs.(2) one derives the Tolman-Oppenheimer-Volkoff (TOV) equa- 
tion that describes nonrotating compact objects in the hydrostatic equilib- 
rium: 



dP{r) _ GM{r)p{r) 

X 




2GM(r) 
c 2 r 



(3) 



Here we denote by M(r) the integral M(r) = Jq p(r')d?r'. Usually M(r) is 
referred to as " the mass within the radius r" . This name is rather mislead- 
ing, as it suggests that only matter energy density is responsible for building 
the neutron star mass M = M{R). It is our prime goal here to calculate the 
gravitational energy contribution to M. In the following we shall be using 
the differential form of the M(r) definition, 

« = VpM. (4) 

dr 

Choosing an appropriate equation of state of dense matter 

P = P(p), (5) 

and the boundary conditions we solve this set of equations for different 
values of central density p c obtaining the whole family of configurations of 
compact objects in hydrostatic equilibrium. 

p(r = 0) = p c and P(r = R) = 0, (6) 

We report here calculations for four equations of state, which we denote 
as A18, Al8+<5v+UIX*, FPR and MS (see Sect. [§. These four EOS have 
been chosen from about twenty EOS we studied solely to make the plots 
clear. The results of our calculations are shown in Fig. ^ 

Numerical values of the neutron stars masses for these equations of states 
are given in the Tabled We show parameters corresponding to the standard 
neutron star mass M = 1A4Mq, and to the maximum mass neutron star, 
which is different for every equation of state. 



2. Gravitational field energy of the compact objects 

Boundary conditions at the surface of the neutron star, r = R, ensure 
that the external metric is just the Schwarzschild metric with the mass M. 
There is no more matter contributions to M at r > R since the matter 
energy density is zero outside the star. Hence the same mass M survives in 
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Fig. 1. Neutron star masses versus central density for four equations of state. 



Equation of state 


Central density [10 15 g/cm 3 ] 


Mass [M ] 


A18 


2.54 


1.44 


A18 


4.23 


1.68 


Al8+5v+UIX* 


1.02 


1.44 


A18+5v+UIX* 


2.74 


2.22 


FPR 


1.36 


1.44 


FPR 


3.38 


1.80 


MS 


1.01 


1.44 


MS 


2.84 


2.04 



Table 1. Computed neutron star masses. For each equation of state the first line 
shows the value of the central density of 1.44 M neutron star, and in the second 
row the central density of the maximum mass neutron star and the value of this 
mass are given. 



the metric functions in the asymptotic region, r — > oo, and M is the total 
gravitational mass of the neutron star (sometimes called the ADM mass) . 

In General Relativity the mass of the star can be split into two parts: the 
matter energy and the gravitational field energy [I], Mc 2 = E Matter +E Field- 
Formally we can write 

Mc 2 = 4vr / Vo° + t%y l2 e Xl2 r 2 dr, (7) 
J o 

where Tq is the energy density component of the matter energy-momentum 
tensor, and i[} is the energy density of the gravitational field pseudotensor. 
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The matter energy can be easily calculated for any given neutron star, 
as all three functions, p(r), A(r), v(r) entering the integral, 

EMatter = 4vrc 2 /* p e v ' 2 e x ' 2 r 2 dr, (8) 
Jo 

are determined by solving Eqs.(2). One should notice that not only the 
integral (jBJ) is uniquely determined, but also matter energy density is com- 
pletely known, p Matter = pe"/ 2 . This energy density includes contributions 
from nuclear and electromagnetic interactions (through p), as well as the 
gravitational interaction (through e u ^ 2 ). It is important to stress that, in the 
Newtonian limit, the generally relativistic matter energy EMatter includes 
twice the (negative) potential energy contribution, as shown by Lynden-Bell 
and Katz 

Using this definition Q we can compute the matter energy for all equa- 
tions of state we study. The results of our calculations are shown in the Fig. 
121 as functions of the total gravitational mass. 




M[M ] 

Fig. 2. Matter energy contribution to the neutron star mass as a function of the 
total gravitational mass. 

One should perhaps stress once again, that for a given equation of state 
and a chosen central density, the matter part of the total mass integral (7) 
is fully determined. For spherically symmetric metric there does not exist 
any ambiguity as far as the integral (8) is concerned, despite the pseudoten- 
sorial nature of the gravitational field energy density if]- The equations (2) 
are fully immune to the particular choice of the gravitational field energy- 
momentum pseudotensor. 
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Having the matter energy EMatter we can easily find the gravitational 
field energy contribution to the neutron star mass, 

E Field = Mc 2 - EMatter = 4vrc 2 J* (e~ x ' 2 - e v ' 2 ) pe x / 2 r 2 dr (9) 

where, formally, we can refer to the function under the integral as the "grav- 
itational field energy density", 

pFiddir) = (e~2 - e5)p(r), (10) 

in units of gjcvr? . Unfortunately, in contrast to the gravitational field energy 
contribution to the neutron star mass, Ep^i^ which is fully determined from 
Eq.Q, this function is not unique. The reason is that t[j is not unique, only 
the integral in Eq.Q has a well defined meaning. 

In the following we calculate the graviational field energy contribution 
to the neutron star mass in terms of often used pseudotensors (which are 
discussed below in some details). From equation ()16|) we have the field 
energy in the Einstein- Tolman prescription, 

(11) 

Using the M0ller pseudotensor ()19|) we find the gravitational field energy in 
the form 

SjfflT _ 12, £ ^ P(r) ^ZZ*ZZI Sir. (12, 

The definitions of the gravitational field energy of the neutron star Eqs. 
©, (|11|) and l|12|l give the same value of the field energy contribution. This 
holds for all equations of state we study as it is shown in Fig. |31 

Now we can compare contributions of the matter energy and the gravita- 
tional field energy to neutron star masses. In Fig. 0]we show both energies 
for all equations of state we use. For low masses of neutron stars the matter 
energy dominates. For "canonical" neutron star mass M = IAAMq the 
gravitational field energy contributes about 20% to the total gravitational 
mass. This contribution increases significantly for maximum mass neutron 
stars, exceeding 40% for the A18+<5v+UIX* EOS. Crossing point of the 
curves in Fig. ^indicates the star for which half of the gravitational mass is 
comprised by energy of its gravitational field. One can see that this occurs 
for three equations of state, albeit for unstable stars. 

For completeness we provide here some formulae for two pseudotensors 
we used in our calculations. We tried a number of gravitational field pseu- 
dotensors, available in the literature. We tried Einstein- Tolman, M0ller, 




Landau-Lifshitz and Weinberg pseudotensors. However, not all of them 
were suitable for our calculations. To obtain the formulae given below we 
calculated all the components in Gallilean coordinates (see Appendix A) and 
employed the GRTensor program to perform the algebra. Eventually, we 
performed numerical integration indicated in Eqs.(8, 9, 11, 12) for Einstein- 
Tolman and M0ller pseudotensors. 



2.1. Einstein- Tolman pseudotensor 

The energy- momentum pseudotensor of Einstein and Tolman |I] is given 
by the formulae: 

^ = 16^- (13) 
where 9" = y/=g% + T»), and 

Hp = [-g (fPgl* - g"><r)\ „ ■ (14) 
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Fig. 4. Comparison of the gravitational field energy and matter energy of neutron 
star 



Here g is the determinant of the metric tensor. This pseudotensor satisfies 
the local conservation laws: 

-± = 0. (15) 
dx u y ' 

For the metric Q we find the energy density for the Einstein- Tolman pseu- 
dotensor to be 

^ _ GM{r) (P(r) + c 2 p(r)) 
*° ~ -2GM(r)+c^r ' 

Some technical details of the calculations are shown in the Appendix B. The 
same result was obtained in Ref. 0. 



2.2. M0ller pseudotensor 
The energy-momentum complex of M0ller ,7 a is given by the expression 

"Z-SRjxT" (17) 
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where = \/—g{t v „ + T%) and the antisymmetric tensor ls 

XV = V=9 (w -9u P ,* ) 9 up 9^. (18) 
For metric (^Q) we have : 

t° = 3P(r). (19) 
Some useful formulae are shown in the Appendix C. 

2.3. Landau pseudotensor 
The Landau pseudotensor reads [S]: 

IT = J^S^, pa (20) 

where S^ up(T = -g{g^ v g p(T - g pp g va ) and IF" = -g{T^ u + t pv ). 
The energy density is then given by 

_ c 2 GM(r) [M(r)-47rr 3 p(r)] 

2nr 3 (2GM(r) - c 2 r) 1 ' 

One can notice that this energy density becomes negative outside the stellar 
radius. This is unphysical and we discard the Landau pseudotensor in the 
following. The fact that the Landau pseudotensor gives the negative en- 
ergy density of the gravitational field of the Schwarzschild metric was first 
recognized by Virbhadra [2- One could probably attribute this behaviour 
to a peculiar dependence on —g, as observed by Cooperstock and Rosen 
jlUj . Some details of derivation of the t[] component of the Landau-Lifshitz 
pseudotensor, are given in Appendix D. 

As far as the Weinberg's pseudotensor is concerned, the relevant com- 
ponents vanish (Appendix E). 



3. The problem of localization of the gravitational field energy 

The fact that all three definitions of gravitational field energy of neutron 
stars give the same value is reassuring. One can be quite confident that 
what we find is really the prediction of the general relativity theory as to 
the gravitational energy field contribution to the neutron star mass. 

There is another important conclusion to be drawn: In all cases we 
consider the gravitational field energy is rather well localized. To show this 
we plot in Figs. |S] and " gravitational field energy density" for all three 
definitions we use. Although these are not to be understood as physically 
meaningful, plots in Figs. |S] and |H1 show they share the same feature: namely 
all densities are localized within the star. 
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The first density corresponds to our definition of the field energy lf§]) 
and is given by Eq.(10). The other two gravitational field energy densities 
correspond to Einstein- Tolman pseudotensor, 

e-t _ e»l*GM(r) (P(r) + c 2 p(r)) 
P Field- -2GM(r)+c 2 r [ ' 



and to the M0ller pseudotensor, 



Moller 
PField 



3e u/2 P(r) (23) 




Fig. 5. The energy density of the neutron star gravitational field according to 
different definitions for the "canonical" mass neutron star 1.44 Mq. 

In Figs. IHllHlwe show results for "canonical" neutron star and for max- 
imum mass neutron stars respectively. 

The fact that the whole gravitational field energy is accounted for within 
the star radius is fully consistent with the findings of Nahmad-Achar and 
Schutz who found that the energy associated with the Schwarzschild 
space-time was independent of radius for radii exceeding the radius of spher- 
ically symmetric mass. We thus arrive at the conclusion that for spherically 
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symmetric compact objects with Einstein- Tolman and M0ller gravitational 
field pseudotensors, the whole gravitational mass is localized within the 
star, including the gravitational field contribution which for neutron stars 
is substantial - reaching 40% of the total mass for maximum mass neutron 
stars. 

In this case the total gravitational mass, M, which is the mass of the 
exterior Schwarzschild metric has all the contributions inside the stellar 
radius, and thus can be used for calculations of such properities as the last 
stable orbit around the neutron star, for sufficiently massive stars. However, 
there is some unusual property of the gravitational field in this approach: 
the field just below the surface of the star has the energy which contributes 
to the stellar mass, whereas above the stellar surface the field does not 
contribute to the neutron star mass: no gravitational field energy is localized 
outside the stellar surface. 
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Fig. 6. The same as in Fig. 5 for maximum mass neutron stars 



4. Equations of state and nuclear models 



In this section equations of state of dense matter, used for calculating 
curves in the above plots, are briefly discussed. 
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The density dependence of the pressure P(p) plays the crucial role in 
calculations of neutron star models. The mass density of matter, p, includes 
kinetic and interaction energy density of the nucleons and rest masses of neu- 
trons and protons pins) = £ (n-N,np) + njymjy + npmp, while the pressure 
of the nuclear fluid is 

The baryon number density np = n^ + np is the sum of the neutron and 
proton number densities. 

We chose four realistic nuclear interaction models which give us the 
energy density as a function of the baryon number density, e{np). These 
are interactions derived by Myers and Swiatecki (MS), the Friedman 
and Pandharipande interaction model XZ i (as parametrized by Ravenhall in 
Ref-HU) (FPR) and two modern potential models, A18 and A18+<5v+UIX* 
|15j . which provide a good fit to the nucleon-nucleon scattering data in the 
Nijmegen data base. 

Myers and Swiatecki proposed a nuclear model based on the Thomas- 
Fermi approximation |161I17| . They extended the model of Seyler and Blan- 
chard |18j by improving four-parameter interaction force between nucleons. 
By fitting two additional parameters they obtained an excellent agreement 
with the ground state energies of about 1600 atomic nuclei. The energy 
density reads 



£ms (n N , n P ) = A\n 3 N + np) - B\ (n 2 N + np 



8 



-2B 2 nNn P + Ci\n 3 N + np 

5\ / 2 5' 

+C 2 I n 3 N n P + n N np — D I 5n^n P - np) , (25) 



where the parameter values are: 



A = 88.247 MeV fm 2 B x = 169.217 MeV fm 3 
B 2 = 286.715 MeV fm 3 C x = 358.664 MeV fm 5 
C 2 = 1167.630 MeV fm 5 D = 50.148 MeV fm 2 . 

D. G. Ravenhall has found a simple formula to reproduce nuclear matter 
calculations of Friedman and Pandharipande ^Hl- The interaction model 
fits the s-, p-, d-, and f-wave phase shifts in the low density limit, and at 
high density the hypernetted chain calculations which makes the effective 
interactions similar to a contact two-body force. Adjustable tensor force 
contribution is also included to obtain good saturation properties. These 
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model forces have been fitted to the nuclear matter calculations at low and 
high densities. The energy density is 



£fpr (nN,np) 



+ B N ) t n + ( + Bp ) T P 



2m n 



+n 2 B 



ai + a 2 e- binB + 



+n B e' b2n B 



(H(< 

05 + aen B + Q - (a 7 + a 8 n B ) 



2mp 
a 3 + a 4 e- fcins 

2 



(26) 



where: 



n B = n N + n P , 



np 



X 

itti 

Bi = (a 9 n B + a w ni) e~ bsnB , n = § (3vr 2 ) t „. 



for i = N,P. The parameter values are: 



res ' 

2 5 
3 (q^2\ 3 3 



01 = 1054.0 MeF/m 3 
a 4 = 2859.0 MeV fm 3 
a 7 = 5.5 MeV 
a w = -59.0 MeV fm 5 
63 = 0.457 fm 3 



a 2 = -1393.0 MeV fm 3 
a 5 = -1.78 MeV 
a 8 = 197.0 MeV fm 3 
h = 0.284 fm 3 



a 3 = -2316.0 MeV fm 3 
o 6 = -52.0 MeV fm 3 
a 9 = 89.8 MeV fm 5 
b 2 = 42.25 fm 6 



Akmal, Pandharipande and Ravenhall used variational chain summation 
methods and created new Argonne v±s two-nucleon interaction which fits all 
of the nucleon-nucleon scattering data in the Nijmegen data base. 

The A18+<k> +UIX* model comprises the boost correction Sv to the two- 
nucleon interaction, which gives the leading relativistic effect of order (v/c) 2 
and the three-nucleon interaction Urbana model IX* which predicts a tran- 
sition in nucleon matter to a phase with neutral pion condensate at the 
density 0.195 fmT 3 for pure neutron matter and 0.32 fvrT 3 for symmetric 
nuclear matter. 

The energy density for both A18 and A18+(5v+UIX* models can be 
given by a common formula, 



£APR {n N ,n P ) 



+ B n )t n + 



+ B p )tp 



,2mjy "J " \2mp 
-An N n P (pi + p 2 n B + p§n B + (p w + pwn B ) e~ p 9 n B^j 

- (n N - n P ) 2 {pviln B + p 7 +psn B + pi 3 e~ p ^ n B^ 
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for A18+<k> +UIX* with rig > p±g or > P20 additionaly: 

-An N n P (n B - P19) (pn + P21 (n B - P19)) e^"™ 9 ) 

- (n N - n P f (n B - p 20 ) (pis + Pi4 (n B ~ P20)) e p (27) 

where: B< = (p 3 n B + p 5 ni) e~ P4 ' riB , n = % (3n 2 )% nj for i = 7V,P. The 
parameters P3 = 89.8 MeU /m 5 , j>4 = 0.457 fm 3 and 775 = — 59.0 MeV fm 5 
are common for these two models. 

For A18 model the parameters are: pi = 297.6 MeV fm 3 , p 2 = -134.6 
MeVfm 6 , p 6 = -15.9 MeV fm 9 , p 7 = 215.0 MeV fm 3 , p 8 = -116.5 
MeVfm 6 , p 9 = 6.42 fm 3 , p w = 51 MeV fm 3 , p u = -35 MeV fm 6 , 
P12 = -0.2 Mey and p 13 = OMeF/m 3 . 

For A18+<k>+UIX* model the parameters are: p\ = 337.2 MeV fm 3 , 
p 2 = -382 MeV fm 6 , p 6 = -19.1 MeV fm 9 , p 7 = 214.6 MeV fm 3 , p 8 = 
-384 MeV fm 6 , p 9 = 6.4 fm 3 , p w = 69 MeV fm 3 , p u = -33 MeV fm 6 , 
P12 = 0.3b MeV, p u = OMeV fm 3 ,p l4 = 0MeVfm 9 ,p 15 = 287 MeV fm 6 , 
pie = -1.54/m 3 ,pi 7 = 175.0 MeV fm 6 , p 18 = -lA5fm 3 ,p l9 = 0.32 fm~ 3 , 
P20 = 0.195 fm~ 3 and p 2 i = MeV fm 9 . 



Appendix A 

Metric in Cartesian coordinates 



First we change the spherical coordinates to the Cartesian ones. Using 
transformation 

r = \/ x 2 + y 2 + z 2 

4> = arctan(y/x) 

9 = &iccos(z/ \j x 2 + y 2 + z 2 ) 

we obtain the metric in Cartesian coordinates, 



ds 2 = e»c 2 dt 2 - eXx [ + y Z + / dx * _ E . , 
x +y z + z z x + y z + z z 

(e A - l)2yz 



e x y 2 + z 2 2 
dy 



+ 



x 2 + y 2 + e x z 2 

x + y 2 + z 2 
( (e x -l)2xy 
\ x + y 2 + z 2 



dz z 



dy 




x + y 2 + z 2 



(A.l) 
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The components of the covariant metric tensor are 



t e v 



V 



o 

e A x 2 +y 2 +z 2 
x 2 +y 2 +z 2 

(-l+e A )it/ 
x 2 +y 2 +z 2 

(— l+e A ) x z 
x 2 +y 2 +z 2 





(-l+e A ) xy 
x 2 +y 2 +z 2 

x 2 +e x y 2 +z 2 
x 2 +y 2 +z 2 

x 2 +y 2 +z 2 



o \ 

x 2 +y 2 +z 2 
— (y z)+e x y z 
x 2 +y 2 +z 2 
x 2 +y 2 +e x z' 2 



(A.2) 



x' 2 +y 2 +z' 2 / 

The determinant of this metric tensor is g = —e x+u . Then we find the 
components of the contravariant metric tensor in Cartesian coordinates to 
be 



( e"" 




V 





x 2 +e x (y 2 +z 2 ) 
e x (x 2 +y 2 +z 2 ) 
(-l+e A ) xy 

e x (x 2 +y 2 +z 2 ) 
(-l+e A ) xz 

e A (x 2 +y' 2 +z 2 ) 





(-l+e A ) x y 
e x (x 2 +y 2 +z 2 ) 
y 2 +e x \x 2 +z 2 ) 
e A (x 2 +y 2 +z 2 ) 
(— l+e A ) yz 
e x (x 2 +y 2 +z 2 ) 



\ 

(-l+e x )xz 

e x (x 2 +y 2 + z 2 ) 

(-l+e x )yz 
e x (x 2 +y 2 +z 2 ) 
e x [x 2 +y 2 )+z 2 
~ e A (x 2 +y 2 +z 2 ) ) 



(A.3) 



Note, that now A and v are functions of Cartesian coordinates x, y and z. 

Appendix B 

Calculation of the Einstein- Tolman pseudotensor 

We are interested in the energy component of the Einstein- Tolman pseu- 
dotensor [3]. Using the equations 1)13(1 and (|14[) one can calculate the non- 
vanishing components of Hq 



H 0l = e (u-X)/2 (3^11) + ^12) + dz{e X g l 3) 

- e ( ^ A)/2 (d x (e X g 21 ) + d y (e x g™) + d z {e x g™) 



tr03 
^0 



>-A)/2 



(dx(e X g 31 ) + d y (e x g 32 ) + d z (e x g 33 



(B.l) 



Inserting the componets of the contravariant metric tensor we have 

2 (e A - l) x 



H 



01 



x 2 + y 2 + z 2 



02 



u 03 
H 



(.» - 1) 



y 



x 2 + y 2 + z 2 
2 (e x - l) z 
x 2 + y 2 + z 2 



(B.2) 



16 PRACA PRINTED ON FEBRUARY 7, 2008 



The 00 componet of 6 pseudotensor is given by 



n — + W + (B.3) 



16vrG 

Taking the partial derivatives of Hq 1 we obtain 



e (A+v)/2 



-l + e A ) P(r) + c 2 (l + e A ) p(r) 



2 < B ' 4 ) 

Notice that e _A = 1 - 2GM(r)/c 2 r, so the 00 component of the Eintein- 
Tolman pseudotensor is 

3 

;l ( -2GM(^r ) 5 (c 4 rKr)+GM(r) (P( r ) - c 2 p(r))) 



90 



(B.5) 



Using definition of ((13)) we obtain the equation (|16() which is the gravita- 
tional field energy density in the Einstein- Tolman prescription. The same 
result was obtained by [Hj. 

_ GM(r) (P(r) + c 2 p(r)) 
to ~ -2GM(r)+c 2 r ' 

Appendix C 

M0ller prescription 

For metric in Cartesian coordinates the x potential is 

xV = -^(^ 9 ^) (CI) 

The nonzero componets of x are 

x 01 = _ e (X + u)/2 " f g U +yg 12 + zg 13 



v' 


V* 2 


+ y 2 + z 2 




v' 




+ y 2 + z 2 




u' 



+ y z + z z v 

x os = _ e (A + „)/2 V ' (.,•,/«' • • zg™) (('.2! 

where prime denotes derivative of the function v with respect to r = 
\/ x 2 + y 2 + z 2 . Inserting the components of the contravariant metric tensor 
and using the equation (|2b|) we have 
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-\ + v 

e 2 x 

,oi _ 



(c 4 (-1 + e A ) + 8 7T G e x P (x 2 + y 2 + z< 



Xo ' c 3 (x 2 + y 2 + z 2 ) 



-X+v 



(c 4 (-1 + e A ) + 8vr G e x P (x 2 + y 
X6~ = " — 



2 i „2 



c 3 (x 2 + y 2 + z 2 ) 



-X+v 

e % z 



. (c 4 (-l + e A )+8^Ge A P(^ 2 + j/ 2 + z^ 

According to the definition (|17[l one finds the 00 component to be 

M ° = ^ (d xX f + d yX ° 2 + a^ 3 ) (C4) 

After taking partial derivatives of potential Xg 7 an d using Eqs. (|2al I2bl l2"cjl 
the M0ller complex reads 



,(A+i/)/2 



-4c 3 GM(r) + 2c 5 r V w KW 

+ e A (-2GAf(r) + c 2 r) (c 4 + 8 GP{r) ir r 2 ) [p{r) + c 2 p{r) 
+ c 6 (5P(r)r-4GAf(r)/o(r)) 

- 4c 4 GP(r) (3M(r) + 2vrr 3 p(r))) (C.5) 
Using the definition of the function A and rearranging terms we obtain 



The M0ller pseudotensor is Mq = y/—g(tQ — T$) and one can easily find 
that in this prescription the gravitational energy density is given by very 
simple formula i(j = 3P(r). 

Appendix D 

Landau- Lifshitz prescription 

For the metric in Cartesian coordinates we have 
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The determinant of the metric tensor is — e I/+A and the contravariant com- 
ponent g 00 is e~ u , so Landau-Lifshitz complex reads 



c 4 



8vrG 



(d x d y e x g 12 + c^eV* + d y d z e x g 23 ) (D.2) 



Inserting the components of contravariant metric tensor, taking partial 
derivatives and using Eqs. (|2al I2bl l2c|) we obtain 



L oo = c 4 (-GM 2 (r) + 2 C 2 7rrV(r)) 
2 vrr 2 (-2 G Af(r) +c 2 r) 2 



For r > Rwe have the Schwarzschild solution. So taking p = we find how 
does the Landau-Lifszyc complex look like in the vacuum outside the star: 
the vacuum energy density is negative. The same result was obtained by 
Virbhadra (see Eq.(16) in [§] with Q = and wiht G = c = 1). Because 
Lq = L 00 goo we finally find the gravitaional field energy density in Landau- 
Lifshitz form to be 

,o = c 2 GM(r) (M(r)-4^r 3 p(r)) 

27rr 3 {2GM(r)-c 2 r) ' 1 ' ' 

Appendix E 

Weinberg pseudotensor 

In the Weinberg prescription |3j the complex of matter and field energy 
density is 

= -^—A^ (E.l) 
167rG 

The tensor A^ ua is given by 

dhP R dhP R dh?v dhP v h^ a h va 

A ,ua = ^ua _ ^ m _ ™rj<° + ™ ^ + 2— - (E.2) 

ox^ ax v ax p oxP ax v ax M 

where = g^ — r]^, and rj^ is Minkowski metric. We look for W 00 . 
Using ho = ri^hpv we find that 

A 00« = K oo _ ^ Oa _ dh^ Oa + dh^ Oa + ^ _ ^ = „ (R3) 

8xq 8xq dx 13 dx 13 8xq Oxq 
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So the Weinberg complex 00 component is equal for any a and (3 as it 
was observed by Xulu [T§] . 
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